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■“^The  potential  difference  can  amount  to  volts  in  the  ionosphere,  and 
kilovolts  in  the  solar  wind,  that  is,  large  compared  with  the  potentials 
typically  generated  by  photoemission  alone.  As  in  the  more  familiar 
case  of  photoelectric  charging,  the  asymmetric  sheath  and  potential 
barrier  produced  by  the  plasma  flow  can  lead  to  erroneous  interpreta- 
tions of  experiments  measurina  space  electric  fields  and  low-energy 
particle  spectra.  In  an  example  of  photoelectric  emission,  a sunlit 
area  on  an  otherwise  dark  surface  becomes  positively  charged  by  the 
emission,  and  is  found  to  acquire  a potential  more  than  twice  the 
emitted  energy  relative  to  the  dark  surface.  This  effect  is  associated 
with  the  physics  of  the ^terminator*^ bounding  the  sunlit  and  dark 
areas . ^ 
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Differential  spacecraft  charging  takes  place  when  the  spacecraft  sur- 
face is  partly  or  entirely  insulating  and  the  charged-particle  fluxes  vary 
from  point  to  point  over  the  surface.  In  the  relatively  familiar  case  of 
photoelectric  emission  from  a sunlit  insulated  area  (cf.  Grard,  1973),  due 
to  electrons  escaping  from  it  the  sunlit  area  tends  to  become  positively 
charged  relative  to  the  surrounding  dark  areas.  Another  mechanism  of  dif- 
ferential charging,  which  is  less  familiar  and  appears  to  have  been 
discussed  little  if  at  all  in  the  literature,  is  that  due  to  the  relative 
motion  between  a nonconducting  spacecraft  and  the  external  plasma  (e.g.,  a 
spacecraft  in  the  ionosphere  or  in  the  solar  wind).  The  fluxes  of  ambient 
ions  and  electrons  on  the  wake  surface  are  not  the  same  as  on  the  front  sur- 
face. For  high  velocities  of  relative  motion  compared  with  the  mean  ion 
thermal  velocity,  whether  this  occurs  in  the  ionosphere  (due  principally  to 
spacecraft  motion)  or  in  the  solar  wind  (due  principally  to  plasma  motion), 
there  is  a significant  differential  in  the  ion  fluxes,  but  a negligible 
differential  for  the  electrons.  Since  the  net  current  density  must  vanish 
locally  at  each  surface  point  in  the  steady  state,  this  plasma-flow  effect 
leads  to  a larger  negative  equilibrium  potential  on  the  wake  surface  than 
on  the  front  surface.  If  there  is  photoemission  as  well  on  the  front  sur- 
face (as  in  the  solar  wind),  this  differential  charging  is  enhanced.  The 
present  report  is  concerned  with  these  effects.  As  shown  below,  differences 
between  the  front  and  wake  surface  potentials  amounting  to  many  kT/e  (where 
T is  the  temperature,  k is  Boltzmann's  constant,  and  e is  the  electron 
charge),  together  with  a potential  barrier  for  electrons,  can  be  generated 
by  the  motion.  The  potential  differential  may  be  expected  to  be  of  the 
order  of  volts  (tens  of  kT)  in  the  ionosphere,  and  of  the  order  of  a kilo- 
volt (100  kT)  in  the  solar  wind.  In  the  ionosphere  and  solar  wind,  this 
differential  can  therefore  be  much  larger  than  that  generated  by  photoemis- 
sion alone.  In  the  magnetosphere,  however,  the  plasma-flow  effect  is  rela- 
tively weak. 


As  of  the  mid-1960's  there  was  already  a considerable  literature  on 
the  subject  of  estimating  spacecraft  potentials,  using  simplified  models 


without  considering  trajectories  and  assuming  perfectly  conducting  space- 
craft (see  reviews  by  Brundin,  1963;  Whipple,  1965;  Samir  and  Will more, 

1966).  These  relatively  crude  models  for  estimating  spacecraft  potentials 
assumed  either  (a)  very  thin  sheaths  such  that  the  "planar  approximation" 
could  be  used,  or  (b)  very  thick  sheaths  (for  small  objects  in  the  iono- 
sphere) but  with  radial  symmetry  so  that  the  simple  Langmuir  theory  could 
be  used.  The  crude  estimates  were  not  concerned  with  differential  charg- 
ing; they  sufficed  for  treating  effects  where  the  detailed  structure  and 
asymmetry  of  the  sheath  (potential  barriers,  for  example)  were  not  con- 
sidered important. 

From  the  mid-1960's  onward,  spacecraft  have  increasingly  sampled  the 
magnetosphere  and  solar  wind,  where  the  spacecraft  conditions  are  altered 
in  several  ways  important  for  differential  charging.  First,  there  is  an 
increase  of  sheath  thickness  to  the  order  of  the  spacecraft  dimensions,  as 
opposed  to  the  thin  sheaths  encountered  in  the  ionosphere;  this  thickness 
is  governed  essentially  by  photoelectrons  and  secondary  electrons  from  the 
surface,  the  space  plasma  contribution  being  typically  relatively  weak. 

A second  circumstance  is  that  many  spacecraft  surfaces  are  partly  or  entirely 
nonconducting  (e.g.,  composed  of  glass-covered  photocells  or  insulating 
thermal  control  blankets),  which  becomes  important  when  the  sheath  is  thick. 

A third  circumstance  is  that  nonuniform  particle  fluxes  occur  over  the 
spacecraft  surfaces,  e.g.,  due  to  photoelectron  and  secondary  emission  and 
to  plasma  flows.  The  combination  of  the  foregoing  circumstances  leads  to 
differential  charging  of  the  spacecraft  surfaces,  which  can  have  deleteri- 
ous effects  as  follows. 

If  the  differential  charging  is  severe  enough,  spacecraft  malfunctions 
can  occur  due  to  electrical  discharges  on  the  insulating  surfaces  (Fredricks 
and  Scarf.  U ^3;  Rosen.  1976) . The  appearance  of  hot  magnetospheric  plasmas 
in  the  kilovolt  temperature  range  impacting  on  the  dark  surface,  combined 
with  photoelectric  emission  on  the  sunlit  surface,  make  possible  intense 
differential  charging  such  that  te^s-of-kilovolts  differentials  can  appear 
(see  spacecraft  charging  symposia  by  Grard.  1973;  Rosen.  1976;  Pike  and 
Lovell.  1977). 
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[ Even  if  the  differential  charging  is  not  so  strong,  say  no  more  than 

\ tens  to  hundreds  of  volts  of  differential,  it  can  interfere  with  measure- 

I ments  of,  say,  weak  ambient  electric  fields  or  low-energy  particle  spectra 

I , (Grard.  1973;  Whipple  and  Parker.  1969) . An  interesting  feature  of  differ- 

ential charging  as  it  affects  low-energy  electron  measurements  is  that  it 
can  create  electron  potential  barriers  which  can  return  emitted  photoelec- 
trons or  secondary  electrons  to  the  surface  and  lead  to  erroneous  interpre- 
tations of  the  data  (Fahleson,  1973) . This  type  of  electron  potential  bar- 
rier is  distinct  from  and  should  not  be  confused  with  the  more  familiar 
space-charge  potential  minimum  of  the  order  of  a volt  which  can  be  produced 
by  ©Hi tted-el ectron  space  charge  and  is  not  due  to  differential  charging. 
The  space-charge  potential  minimum  has  been  studied  theoretically  by  Soop 
(1972,  1973),  Schroder  (1973),  Parker  (1976b),  Whipple  (1976a),  and  Roth- 
well  et  al . (1977);  it  can,  however,  be  less  important  than  the  barriers 
produced  by  differential  charging  effects.  Discussing  a well -documented 
experiment  on  the  ATS-6  geosynchronous  satellite,  Whipple  (1976b)  infers 
that  photoelectrons  and  secondary  electrons  from  the  spacecraft  surfaces 
are  reflected  from  a potential  barrier  which  is  much  too  large  to  be  due 
to  space  charge  but  must  be  associated  with  some  kind  of  differential 
charging  (Whipple,  1976ab).  A similar  potential  barrier  due  to  differen- 
tial charge  is  that  produced  artificially  by  an  attractive  electron  trap 
mounted  on  a repulsive  spacecraft,  which  can  cause  secondary-emission 
currents  to  be  incorrectly  interpreted  as  plasma  currents  (Whipple  and 
Parker,  1969) . In  the  present  report,  a potential  differential  and  a 
potential  barrier  are  shown  to  be  producible  by  a plasma  flow.  This  can 
lead  to  difficulties  of  interpretation  of  solar  wind  measurements  such  as 
those  of  Rosenbauer  (1973),  and  may  be  responsible  for  singularities 
observed  by  photoelectron  detectors  in  the  ionosphere  (W.  K.  Peterson, 
private  communication,  1977). 

The  procedure  for  theoretically  predicting  sheath  asymmetries  and 
, potential  barriers  is  generally  complicated  in  that  it  requires  particle 

trajectory  calculations  as  well  as  a three-dimensional  sheath  description 
for  a realistic  treatment  (Parker.  1970,  1973.  1976a.  1977;  Parker  and 


Whipple.  1967,  1970;  and  appendix  of  this  report). 


In  this  report  we  present  results  of  sample  calculations  of  differen- 
tial charging  due  to  both  plasma  flow  and  photoemission,  primarily  address- 
ing the  asymmetric  sheath  and  potential  barrier  produced  by  the  plasma  flow. 
These  results  may  be  considered  preliminary,  because  the  photoemission  is 
considered  separately  rather  than  simultaneously.  However,  the  differential 
charging  is  enhanced  by  photoelectric  emission  on  the  front  surface.  Space 
charge  is  neglected  compared  with  surface  charges  as  sources  of  the  field; 
this  neglect  has  been  shown  to  be  justified  (Soop,  1973) . Hence,  while  the 
predictions  may  not  be  quantitative  for  an  actual  spacecraft,  they  are  con- 
servative and  indicate  what  may  be  expected:  (a)  in  the  ionosphere  for 
small  insulated  objects,  small  meteoroids,  or  small  parts  of  a spacecraft 
(e.g.,  a painted  antenna)  located  within  the  wake  region  of  a moving  space- 
craft, and  (b)  in  the  solar  wind  for  an  entire  spacecraft. 

In  the  example  of  photoelectric  emission,  a sunlit  area  on  an  other- 
wise dark  surface  becomes  positively  charged  by  emission,  and  is  found  to 
acquire  a potential  more  than  twice  the  emitted  energy.  This  effect  is 
associated  with  the  physics  of  the  "terminator"  bounding  the  sunlit  and  dark 
areas . 


In  Section  2 we  indicate  briefly  the  nature  of  the  numerical  methods 
used,  which  are  presented  in  more  detail  in  Appendix  A. 

Results  are  discussed  in  Section  3. 

A listing  of  the  computer  program  used  is  given  in  Appendix  B. 
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2.  NUMERICAL  APPROACH 


The  model  chosen  to  represent  a spacecraft  is  a truncated  cylinder  of 
approximately  equal  height  and  diameter,  as  shown  in  Fig.  1 where  it  is 
called  (for  brevity)  3i,  "pillbox."  It  is  assumed  that  the  spacecraft  sur- 
face charge  and  the  electric  field  around  the  spacecraft  are  axially  sym- 
metric, and  that  the  electric  field  in  space  is  given  by  the  solution  of 
Laplace's  equation  on  a grid  of  points  in  r-z  space  (including  the  pillbox- 
shaped spacecraft).  This  means  that  when  the  potential  distribution  on  the 
grid  points  has  been  computed  (cf.  Parker  and  Whipple,  1970),  the  potential 
and  electric  field  at  any  point  in  space  may  be  obtained  by  interpolation. 
This  allows  ion  and  electron  trajectories  to  be  computed  within  the  region 
of  space  spanned  by  the  grid.  The  outer  boundary  of  the  grid  represents 
"infinity."  The  method  used  for  computing  the  field  allows  the  use  of  vari- 
able grid  intervals  (see  Parker,  1976a.  1977;  Parker  and  Whipple,  1970;  and 
the  appendix  of  this  report). 

The  normal  component  of  the  ion  or  electron  flux  may  be  computed  at 
any  point  on  the  pillbox  surface  by  evaluating  a triple  integral  in  veloc- 
ity space,  and  following  trajectories  backward  in  time  to  determine  their 
origin  and  therefore  the  value  of  the  integrand.  This  represents  the 
"inside-out"  method  originated  by  Parker  (1964) . Details  are  given  by 
Parker  and  Whipple  (1970).  and  more  generally  by  Parker  (1976a.  1977)  and 
in  the  appendix  of  this  report.  An  alternative  approach  is  the  outside-in 
method  (see  appendix).  For  the  calculations  to  be  discussed,  the  conditions 
at  the  spacecraft  surface  are  represented  by  a discrete  set  of  grid  points 
and  associated  surface  areas,  as  illustrated  (by  12  points)  in  Fig.  2. 

For  the  present  purposes,  the  differential  potential  and  charge  dis- 
tributions on  the  spacecraft  surface  may  be  determined  by  two  different 
approaches.  In  one  approach  we  may  determine  the  potential  at  each  local 
grid  point  by  a relaxation  process  until  the  net  current  density  is  zero. 
This  involves  "cutting-and- trying",  whereby  the  surface  potentials  (12 
values  as  illustrated  in  Fig.  2)  are  first  given  assumed  values,  and  later 
successively  corrected  in  accord  with  the  signs  and  magnitudes  of  the 
resulting  set  of  net  current  densities.  The  surface  potentials  represent 
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the  "inner"  boundary  condition.  In  the  other  approach,  a time-dependent 
method  is  used  whereby  the  surface  charges  at  the  surface  grid  points  are 
updated  after  each  time-step,  during  which  the  net  current  densities  are 
calculated.  The  next  potential  distribution  is  computed  using  the  new 
surface  charges  as  sources.  In  the  computations  to  be  described  next,  a 
grid  of  about  400  points  was  used,  of  which  24  were  used  to  represent  the 
spacecraft  surface  points.  Further  details  are  given  in  the  appendix. 
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3.  SAMPLE  SOLUTIONS 
A.  Plasma  Flow 

In  one  of  the  problems  treated  here,  we  assume  the  nonconducting  space- 
craft to  be  in  a flowing  plasma,  with  the  plasma  flow  along  the  axis,  from 
the  bottom  (front  region)  toward  the  top  (wake  region)  in  Fig.  1.  The 
plasma  is  taken  to  be  ionized  hydrogen  and  is  assumed  to  have  a velocity  of 
flow  four  times  larger  than  the  most  probable  ion  thermal  velocity  (ion 
"Mach  number"  = 4).  Since  the  unperturbed  ion  flux  to  the  wake  surface  is 
about  9 orders  of  magnitude  smaller  than  the  corresponding  ion  flux  to  the 
front  surface,  and  since  the  electron  fluxes  are  about  the  same  to  the 
front  and  wake  surfaces,  there  will  be  a significant  differential  between 
the  equilibrium  potentials  at  the  front  and  wake  surfaces  (see  below).  The 
effects  of  photoemission  for  the  pillbox  geometry  are  treated  in  the  next 
section. 

Figure  3 shows  equipotential  contours  around  the  spacecraft,  obtained 
by  numerical  solution,  labeled  by  encircled  numbers  representing  dimension- 
less values  of  the  potential  (in  units  of  kT/e,  where  T is  the  plasma  tem- 
perature). The  errors  in  the  solution  shown  are  estimated  to  be  under  10 
percent.  These  potentials  are  obtained  from  Laplace's  equation,  where  the 
surface  potentials  are  obtained  by  the  relaxation  method  discussed  in  the 
appendix,  under  the  requirement  of  zero  net  current  density  at  all  surface 
points. 

There  are  three  regions  of  characteristic  behavior  of  the  potential, 
the  "top"  or  "wake",  the  "side",  and  the  "bottom"  or  "front".  In  the  top 
region,  the  potentials  are  of  the  order  of  -lOkT/e.  This  large  negative 
value  (about  -17kT/e  at  the  surface)  is  associated  with  the  reduction  in 
ion  flux  due  to  the  flow.  In  the  side  region,  the  potentials  are  of  the 
order  of  -3kT/e*,  this  is  essentially  the  order  of  the  equilibrium  potential 
when  there  is  no  flow  (--In/m^./mg  kT/e).  In  the  bottom  region,  the  poten- 
tials are  of  the  order  of  -kT/e,  i.e.,  less  negative  than  the  side,  because 
of  the  enhancement  of  the  ion  flux  due  to  the  flow.  (Adding  photoemission 

j here  would  make  the  front  potential  still  less  negative.)  The  surface 
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points  are  thus  not  equipotential . Note  that  there  is  a saddle-point  in 
the  front  region,  that  is,  a potential  barrier  for  electrons.  The  barrier 
height  is  between  10  percent  and  20  percent  higher  than  the  potentials  at 
the  nearest  (front)  surface  points.  This  feature  is  caused  by  the  interac- 
tion between  the  relatively-large-magnitude  wake  potentials  and  the  rela- 
tively-! ow-magnitude  front  potentials.  The  dashed  part  of  the  contour 
labeled  "-3.0"  near  the  side  surface  indicates  that  there  is  more  compli- 
cated fine  structure  (variations  of  potential  along  the  side  surface)  than 
is  shown  in  the  figure.  The  potentials  along  the  top  surface  fall  off  to 
the  right.  The  potentials  along  the  bottom  surface  first  fall  with  radius, 
then  rise  sharply  as  the  corner  is  approached.  This  may  be  a "corner 
effect." 

On  the  basis  of  these  results,  one  would  expect  a relatively  small 
body  in  the  ionosphere,  such  as  a thin  antenna  or  boom  painted  with  noncon- 
ducting paint,  or  a painted  or  insulated  object  in  the  very  near  wake  of  a 
spacecraft  (or  the  spacecraft  surface  itself  if  it  is  a dielectric)  to 
become  highly  negatively  charged,  to  potentials  of  the  order  of  volts  in 
the  ionosphere. 

In  the  solar  wind,  the  above  calculation  could  apply  to  an  entire 
spacecraft  since  the  Debye  length  is  large.  However,  the  ion  Mach  numbers 
are  of  the  order  of  10  rather  than  4,  which  would  lead  to  negative  dimen- 
sionless wake  potentials  an  order  of  magnitude  larger  than  the  wake  poten- 
tials shown  in  Fig.  3 (It  is  shown  in  the  appendix  that  the  potential 
difference  in  units  of  kT/e  generated  by  the  flow  should  be  of  the  order 
of  the  square  of  the  ion  Mach  number  of  the  flow.)  This  means  that  for 
T=10ev  in  the  solar  wind,  one  may  have  kilovolt  potential  differences 
between  the  wake  and  front  surfaces.  The  electric  fields  due  to  this  dif- 
ferential charging  may  significantly  disturb  measurements  of  low-energy 
plasma  electrons,  for  example  on  the  HELIOS  spacecraft  (Rosenbauer.  1973) . 
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B.  Photoemission 


In  the  second  problem  treated  here,  the  bottom  surface  of  the  pillbox 
is  assumed  to  be  illuminated  by  the  sun  (along  the  axial  direction)  and  to 
emit  photoelectrons,  while  the  sides  and  top  of  the  pillbox  are  dark  and 
nonemitting.  (The  axial  direction  of  illumination  is  appropriate  for  main- 
taining axial  symmetry.  In  the  solar  wind  the  plasma  flow  and  solar  illumi- 
nation are  in  the  same  direction.)  The  ambient  plasma  contributions  are 
not  considered  simultaneously  with  the  photoemission.  The  photoelectrons 
are  assumed  to  be  emitted  isotropically  and  monoenergetically,  with  1 ev  of  ’ 
kinetic  energy.  All  points  and  their  associated  areas  on  the  bottom  surface 
are  emitting  except  for  the  corner  point,  for  instance.  Nos.  10,  11,  and  12 
in  Fig.  2.  In  the  actual  problem  the  "terminator"  was  put  at  R®0.95R^. 

(The  terminator  is  not  put  exactly  at  the  corner  R=Rq»  for  numerical  reasons.) 
The  time-dependent  method  is  used,  together  with  the  outside-in  method  dis- 
cussed in  the  appendix.  At  zero  time,  there  is  no  charge  on  any  surface. 

As  time  increases  from  zero,  emitted  photoelectrons  from  the  surface  at 
first  escape  to  infinity,  leaving  behind  positive  charge  and  causing  the 
bottom  surface  to  acquire  a positive  potential.  As  this  potential  builds 
up  to  a value  of  the  order  of  a volt  (the  ejection  energy),  the  photoelec- 
trons no  longer  all  escape  to  infinity,  but  begin  to  return  to  the  space- 
craft. (It  is  of  interest  that  the  potential  is  found  to  rise  significantly 
above  one  volt  - see  below.) 

Figure  4 shows  the  time- behavior  of  the  potential  of  the  center  point 

on  the  bottom  surface  (for  instance.  No.  12  in  Fig.  2).  The  potential  in 

volts  is  plotted  against  dimensionless  time,  where  the  scale  time  t^  in 

seconds  may  be  written  as  1.11  Vq/(RqJq),  where  is  the  kinetic  energy 

of  emission  in  volts,  R.  is  the  scale  dimension  of  the  emitting  area  in  cm, 

0 2 

and  is  the  emitted  current  density  in  pi  coamp  per  cm  . Thus,  for  a 
spacecraft  radius  Rq*50  cm,  photoemission  current  density  Jq»1000  picoamp/ 
cm^,  and  emission  energy  V^=l  volt,  the  scale  time  is  t^^Z.ZZ  x 10"^  sec. 

(The  ordinate  in  Fig.  4 scales  with  the  ejection  energy.)  The  potential 
along  the  bottom  surface  is  not  uniform,  but  is  maximum  in  the  center  and 
falls  off  with  radius  (see  Fig.  5,  which  indicates  an  approximate  drop  of 
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FIG.  4.  ROTENTIAL  OF  CENTER  POINT  VERSUS  DIMENSIONLESS  TIME 
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IG.  5 . DIFFERENTIAL  CHARGING  OF  NONCONDUCTING  SPACECRAFT  BY  PHOTOEMISSION 

- EQUIPOTENTIAL  CONTOURS  IN  UNITS  OF  EMISSION  KINETIC  ENERGY  - time»2t 


30  percent  over  the  illuminated  area).  The  central  potential  rises  to  an 
approximate  plateau  of  slightly  over  2 volts  in  about  one  characteristic 
time  t^,  and,  after  an  interval  of  approximately  constant  behavior  between 
tjj  and  3t^,  continues  to  increase  but  much  more  slowly  than  the  initial  rate. 

It  is  a curious  fact  that  the  potential  should  rise  to  more  than  2 
times  the  ejected  energy  potential -equivalent,  rather  than  to  exactly  the 
ejected  energy  as  may  be  expected  purely  on  the  basis  of  conservation  of 
energy  without  consideration  of  the  nature  of  the  trajectories,  surface,  and 
potential  distribution.  The  computer  results  show  a transfer  of  electrons 
from  the  illuminated  areas  to  the  small  dark  area  associated  with  the  corner 
point;  that  is,  emitted  electrons  are  pulled  back,  but  cross  the  local  ter- 
minator at  R=0.95R^  and  hit  the  dark  corner  area  rather  than  return  to  other 
points  on  the  emitting  area.  The  charge  on  this  corner  area  continues  to 
increase  negatively,  while  the  illuminated- surf ace  charge  increases  posi- 
tively. No  charges  are  deposited  on  the  top  and  side  dark  surfaces.  A 
similar  build-up  to  more  than  the  ejection  energy  has  been  observed  in  com- 
puter experiments  performed  by  De  and  Criswell  (1977)  and  by  Pelizzari  and 
Criswell  ( 1 977 ) in  studies  of  the  photoelectric  charging  of  locally-sunlit 
areas  in  the  dark  lunar  terminator  region.  A possible  explanation  of  this 
“excess"  charging  phenomenon  is  proposed  here  as  follows. 

The  effect  is  appropriate  to  the  problem  of  electron  emission  from  a 
restricted  area  of  a nonconducting  surface,  with  no  compensating  electron 
flux  from  an  external  plasma.  It  also  depends  on  the  returning  electrons 
"sticking"  where  they  hit.  After  the  initial  potential  buildup,  despite 
the  deposition  of  negative  charge  on  the  dark  side  of  the  terminator,  the 
surface  potential  falls  off  monotonically  from  the  central  value  but  remains 
positive  as  one  goes  into  the  dark  region.  That  is,  the  surface  gradient 
(tangential  electric  field)  remains  finite  and  continuous  across  the  ter- 
minator. (If  the  emitting  area  were  a conductor,  the  surface  gradient 
would  be  discontinuous  and  singular  across  the  terminator.)  Thus,  there  is 
a finite  interval  straddling  the  terminator,  such  that  within  this  interval 
electrons  can  cross  the  terminator  from  a sunlit  point  (moving  "uphill")  to 
a dark  point  where  they  are  held  fast,  without  using  up  their  kinetic  energy. 
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A finite  rate  of  transfer  across  the  terminator  from  the  sunlit  area  to 
the  dark  area  can  thus  occur  as  long  as  the  surface  potential  gradient  is 
finite  at  the  position  of  the  terminator,  regardless  of  how  high  the  cen- 
tral potential  of  the  sunlit  area  becomes.  The  transfer  rate  should 
approach  zero  as  the  gradient  at  the  terminator  approaches  infinity. 

Whether  this  process  is  self-limiting,  that  is,  whether  the  gradient  at 
the  terminator  becomes  infinite  within  a finite  time,  is  presently  unknown. 
The  key  point  is  probably  that  the  sunlit  area  cannot  be  strictly  an  equi- 
potential  surface.  It  may  be  approximately  so  over  most  of  its  area,  due 
to  electron  transport  tending  to  maintain  equi potentiality  (De  and  Criswell. 
1977;  Pelizzari  and  Criswell.  1977) . but  the  potential  should  fall  off  in 
the  vicinity  of  the  terminator. 

Figure  5 shows  a few  equi potential  contours  around  the  spacecraft. 

(Only  half  of  the  spacecraft  is  shown,  since  it  is  symmetric  about  the  axis.) 
The  potential  contours  are  taken  from  the  solution  of  the  foregoing  problem, 
at  the  time  t*2t^.  At  this  time  the  bottom-surface  potential  is  approxi- 
mately at  its  plateau  value  (Fig.  4)  of  about  2 volts,  and  is  undergoing  its 
smallest  rate  of  change.  The  source  of  the  field  is  a nonuniform  disk  of 
charge  at  the  bottom  surface,  with  positive  charge  on  the  left  side  of  the 
terminator,  and  negative  charge  on  the  right.  The  equipotentials  are  sym- 
metric about  the  plane  of  the  disk  since  the  spacecraft  dielectric  constant 
was  assumed  to  be  equal  to  the  free-space  value  (no  polarization  effects). 


APPENDIX  A. 


COMPUTER  APPROACHES 

The  computer  program  listed  in  Appendix  B was  developed  for  the  study 
of  the  sheath  about  a pillbox-shaped  spacecraft.  The  sources  of  the  sheath 
electric  field  include  spacecraft  surface  potentials  and  space  charge  due 
to  the  ambient  plasma.  The  program  has  options  to  include  the  effects  of 
insulating  as  well  as  conducting  spacecraft  surfaces,  and  of  electron  emis- 
sion due  uO  photoelectric  or  (with  minor  modifications)  secondary  processes. 
The  program  can  solve  a coupled  Poisson-Vlasov  system  of  nonlinear  partial- 
differential-  integral  equations,  and  uses  a special  iteration  algorithm  to 
obtain  self-consistency  between  the  Poisson  and  Vlasov  solutions.  This 
yields  distributions  of  electric  potential  and  ion  and  electron  density. 

The  "inside-out"  method,  which  follows  ion  and  electron  trajectories  back- 
ward to  their  origin  at  the  body  surface  or  in  the  undisturbed  plasma,  may 
be  used  as  one  option  to  compute  the  necessary  moment-integrals  for  particle 
density  and  flux  at  arbitrarily  chosen  points.  This  is  useful  for  contri- 
butions from  the  plasma.  An  "outside- in"  approach  is  a useful  alternative 
option  for  contributions  from  surface  emission.  A grid  is  used  to  define 
the  spatial  distributions  of  potential  and  density  in  the  space  about  the 
spacecraft.  The  approach  is  applicable  to  a larger  range  of  the  parameters 
than  other  available  approaches. 

Figure  A1  illustrates  the  nature  of  the  grid  representation  used  for  a 
satellite  having  the  form  of  a finite- length  cylinder.  The  geometry  is 
axi ally-symmetric,  with  the  axis  shown  as  the  vertical  dotted  boundary  line 
on  the  left.  The  boundary  condition  representing  the  condition  at  infinity 
is  applied  to  the  other  boundary  lines  of  the  grid.  There  is  an  inner  boun- 
dary representing  the  satellite  surface,  on  the  grid  points  of  which  the 
surface  potentials  ((p^,  etc.)  are  defined.  The  grid  points  are  at  the 
intersections  of  the  grid  lines  at  constant  r and  constant  z.  Associated 
with  each  grid  point  is  a volume  of  revolution  in  the  shape  of  a torus  of 
rectangular  cross-section  (shown  as  shaded  boxes  surrounding  grid  points). 
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The  particle  fluxes  calculated  at  the  satellite-surface  grid  points  I 

I determine  the  equilibrium  surface-potential  distribution  etc.),  | 

with  current  density  balance  at  individual  grid  points  in  the  case  of  I 

? insulating  sections,  and  current  balance  over  the  conducting  areas  in  the 

I ' case  of  conducting  sections.  The  flux  calculation  is  based  on  the  following. 

j . ...  . . 

t Assuming  a fixed  electric  field  configuration  (potential -values  given 

! . ] 
at  all  grid  points),  the  flux  may  be  written  as  a triple  integral  over  ] 

! velocity  space,  of  the  form  | 


j(r)  = I I j f(r,v)  v^  d^ 


(A-1) 


where  v is  the  vector  velocity  of  a particle  passing  through  the  point  r, 
and  where  v^  is  the  component  of  v normal  to  the  surface  at  r.  The  dis- 
tribution function  f(r,v)  is  the  density  of  points  in  six-dimensional 
phase  space.  In  the  absence  of  collisions  and  time-variations,  f depends 
only  on  the  constants  of  motion  and  is  constant  along  any  orbit. 

The  region  of  interest  may  be  considered  to  be  enclosed  by  a composite 
source  surface  on  all  points  of  which  f is  assumed  to  be  known.  There  is 
an  external  boundary  surface  at  infinity  where  f has  the  unperturbed  ambient 
value  f„(v_^),  and  internal  source  surfaces  on  the  spacecraft  where  f has  the 
value  in  accord  with  the  emission  from  these  surfaces.  Separate 

contributions  to  the  flux  j from  infinity  and  from  the  spacecraft  surfaces 
can  be  written  as  independent  integrals  of  the  form  of  Eq.  (A-1),  where  each 
is  comprised  only  of  contributions  from  its  associated  source  surface.  In 
order  to  determine  whether  a specified  velocity  v (at  f)  connects  with  infin- 
ity or  with  a point  elsewhere  on  the  spacecraft  surface,  the  orbit  is  fol- 
lowed backward  in  time  to  its  source  (all  orbits  being  dynamically  reversible). 
This  is  generally  an  appropriate  task  for  a computer  and  is  the  heart  of  the 
"inside-out"  method. 

Since  symmetric  velocity  distributions  are  of  interest  the  polar  form 
of  velocity  space  is  used  in  Eq.  (A-1).  Thus,  the  flux  at  a surface  point  r 
may  be  written  as 
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(A-2) 


r 


j(r)  = 


f V dv  cosa  dQ  (A-2; 

heml:^.iere 

where  the  velocity- space  volume  element  has  been  expressed  in  terms  of  a 
local  velocity-coordinate  system  by 


2 

d*^  = V dv  d£2,  dSJ  » sina  da  dB 


(A-3) 


1 


Here,  v,  a,  and  B denote  the  magnitude,  polar  angle,  and  azimuthal  angle, 
respectively,  of  the  velocity- vector  v,  where  the  definitions  of  the  angle 
variables  a and  B are  illustrated  in  Fig.  A2.  The  flux  j is  the  component 
of  the  flux  vector  j in  the  direction  of  the  chosen  axis,  e.g.,  the  normal 
to  the  surface  at  the  point  r.  In  Eq.  (A-2)  the  subscript  "hemisphere" 
denotes  that  the  angular  integration  is  over  the  hemisphere  of  outgoing 
directions  (2^  steradians). 

Contributions  from  Infinity 

Assuming  at  infinity  a Maxwellian-with-drift  velocity  distribution, 
with  temperature  T,  particle  mass  m,  particle  density  n^,  and  drift  Mach 
number  M,  the  flux  contribution  from  ambient  particles  may  be  written: 


1/2 


J.  • "o  (s)  L 


'Max((^,o) 


f g-U  & cosa  do 

ir 


where 


» E + COS9. 


(A-5) 


Here,  « is  a "cut-off"  function  which  is  unity  if  the  orbit  connects  with 
infinity,  and  is  zero  otherwise;  E and  4)  denote  local  dimensionless  total 
energy  and  electric  potential,  normalized  by  kT  and  kT/e,  respectively; 

9^  is  the  angle  between  v^  and  the  drift  direction  at  infinity,  where  v^ 
is  the  velocity  at  infinity  of  the  orbit  characterized  locally  by  E,  a, 
and  B. 
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Contributions  from  Surface  (Inside-Out  Method) 

A.  Isotropic  Emission 

Assuming  that  at  a surface  point  r there  is  an  isotropic  Maxwellian 
distribution  of  emitted  kinetic  energies,  corresponding  to  surface  parti- 
cle density  n^^  and  surface  temperature  T^,  the  flux  contribution  from  the 
s-th  surface  on  the  spacecraft  may  be  written 


cosa  dfi 


'^s'min 


dE  (E-t)  I I (A-6) 


where 


^so^'s^ 


^ '^so^’^s^ 


kT^Cr^) 


(A-7) 


represents  the  unperturbed  flux.  Here,  5^  is  the  cut-off  function  for  the 
s-th  surface,  i.e.,  unity  if  the  orbit  connects  with  the  s-th  surface,  and 
zero  otherwise;  E and  (|)  denote  local  dimensionless  total  energy  and  elec- 
tric potential,  normalized  by  kT^  and  kT^/e,  respectively. 

B.  Cosine  Emission 

If  the  emission  is  Maxwellian  in  energy  but  proportional  to  the  cosine 
of  the  angle  made  with  the  surface  normal  rather  than  isotropic,  then  Eq. 
(A-6)  is  replaced  by 


dE  (E-*)  f f \3cos.^^dii 

J J 


where  the  function  is  still  given  by  Eq.  (A-7),  but  where  (4/3)j, 
represents  the  new  unperturbed  flux. 
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Contributions  from  Surface  (Outside- In  Method) 


One  can  also  use  an  "outside-in"  approach,  i.e.,  following  trajectories 
forward  in  time,  to  compute  fluxes  or  densities.  This  approach  is  of  course 
equivalent  to  and  would  provide  the  same  results  as  the  inside-out  method, 
but  may  be  computationally  more  efficient  in  certain  problems.  In  particu- 
lar, it  should  be  more  efficient  for  calculating  flux  contributions  from 
the  surface.  For  this  approach,  let  dA  represent  an  element  of  surface 
receiving  contributions  from  another  element  of  surface  dA’.  Moreover,  let 


dpAAi  denote  the  fraction  of  emitted  trajectories  from  dA*  which  reach  dA. 


Morover,  let  j^  denote  the  flux  received  at  A,  while  j^,  denotes  the  flux 
emitted  at  A*.  Then  by  conservation  of  particles  we  may  write 


dj^  • dA  = dF 


AA' 


dA' 


(A-9) 


Thus,  the  flux  received  at  point  A due  to  emission  from  all  emitting  points 
A'  is  given  by  the  integral 


I j 

Jail  A' 


dF 


AA' 


A' 


dA 


dA' 


(A-10) 


where  dF/dA  is  determined  by  trajectory  calculations.  The  foregoing  inte- 
grals may  be  implemented  numerically  as  shown  below. 


Assuming  that  the  ambient  plasma  flows  along  the  z-direction  with  ion 
Mach  number  M,  the  dimensionless  velocity-distribution  function  at  infinity 
may  be  written: 


1 


''o  ° -in 

TT 


(A-ll) 


where  (repeating  Eq.  (A-5)) 


U = E + M^  - ZME^''^  cose 


(A-5) 
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The  flux  integral,  namely 


j(r) 


f(f,v)  v^d^v 


(A-1) 


where  f is  the  distribution  function  and  is  the  component  of  the  particle 
velocity  normal  to  the  surface,  may  be  approximated  for  numerical  evaluation 
by  the  inside-out  method  using  a quadrature  sum  as  follows: 


J * 


I J K 

IZIA 

i j k 


(E 


ijk  '^ijk 


k-i*) 


irr  jk^ 


(A-12) 


where  U is  defined  by  Eq.  (A-5),  (fr  is  the  dimensionless  local  potential, 
and  the  3 indices  refer  to  discrete  values  of  the  3 components  of  velocity. 
The  discrete  values  are  chosen  in  accord  with  a Gaussian  quadrature  scheme, 
and  the  coefficients  A^jj^  are  proportional  to  the  associated  weights  and 
other  factors. 

The  flux  integral  for  the  outside-in  method,  namely,  Eq.  (A-10),  may 
be  evaluated  by  a four- fold  quadrature  sum  as  follows: 


J * 


I J K 

I ^ I Z B--U  • J 

h t ^ ijka 
1 J k a 


A 


4-  . A. 


(A-13) 


Here,  in  addition  to  the  sum  over  the  3 velocity  indices,  we  indicate  a 
summation  over  contributing  finite  areas  aJ  . The  area  A denotes  the  finite 

a 

area  at  the  point  of  interest;  F,  denotes  the  fraction  of  emitted  trajector- 

ies  (assumed  all  emitted  from  the  center  of  A,  ) which  intersect  A;  j, 

1 ® “ 
denotes  the  emitted  flux  at  ; and  is  the  appropriate  coefficient. 

For  the  purpose  of  computing  surface-potential  distributions  over  a 
spacecraft  surface  with  a variation  of  surface  properties,  we  associate 
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separate  differential  area  segments  with  the  surface  grid  points  defining 
the  surface,  as  illustrated  in  Fig.  2.  The  areas  associated  with  the 
surface  grid  points  in  Fig.  2 are  labeled  A2.  etc.  These  may  repre- 
sent areas  with  different  conductivities  and  different  emission  properties, 
either  intrinsic  or  depending  on  their  geometric  position  and  orientation 
(as  in  photoemission). 

For  a given  distribution  of  surface  properties,  the  eoui librium  poten- 
tial distribution  may  be  determined  by  relaxation  as  follows.  A distribu- 
tion of  surface  potentials  (at  the  grid  points)  is  initially  assumed.  This 
leads  by  trajectory  computations  to  a distribution  of  fluxes  of  ions  and 
electrons  at  the  surfaces.  The  ion  and  electron  fluxes  can  be  determined 
numerically  using  surrmations  represented  by  Eq.  (A-12)  or  (A-13).  On  each 
conducting  area  (or  collection  of  connected  conducting  areas)  the  potential 
must  be  adjusted  so  that  the  net  current  to  the  area  is  zero.  At  nonconduct- 
ing points  the  net  current  density  must  be  zero.  These  adjustments  can  be 
accomplished  by  an  iterative  relaxation  procedure.  Thus,  the  initial  guesses 
for  the  surface  potentials  are  modified  according  to  the  sign  and  magnitude 
of  the  net  current  or  current  density.  The  essence  of  the  relaxation  procedure, 
for  a totally  nonconducting  spacecraft,  is  the  iterative  algorithm: 


+ a-CD' 


(A-14) 


where  the  superscript  N denotes  the  N-th  iteration,  and  where  is  the 
dimensionless  potential  at  the  k-th  surface  point,  CDj^  is  the  net  current 
density  at  that  point,  and  a is  the  relaxation  parameter. 

For  the  time-dependent  approach,  the  charge  associated  with  the  finite 
area  A is  given  by  the  time-integral : 

- I j . A . dt  (A-15) 

This  charge  is  used  as  the  source  term  at  the  center  of  the  associated  vol- 
ume of  revolution  indicated  in  Fig.  A1 . In  a step-by-step  procedure,  during 
one  half  of  each  cycle,  the  fluxes  are  evaluated  in  accord  with  Eq.  (A^12) 


25 


or  (A-13);  in  the  other  half  of  the  cycle  the  charges  are  updated  in  accord 
with  the  discretized  form  of  Eq.  (A-15). 


Charging  of  the  Wake  Surface  in  Plasma  Flow 

One  may  estimate  the  potential  acquired  by  the  wake  surface  of  a non- 
conducting spacecraft  in  a plasma  flow,  as  follows.  Assume  that  the  dimension- 
less potential  ♦ is  negative,  and  that  the  electron  and  ion  fluxes  to  the  wake 
surface  are  given,  respectively,  by 

jg  » n^  /kT/Zimig  exp{+^)  {A-16) 

and 

- n^  AT/Zm.  {exp(-M^)/2M2}  (A-17) 

where  M is  the  ion  Mach  number.  For  large  values  of  M,  the  bracketed  expression 
in  Eq.  (A-17)  is  the  asymptotic  form  of  {exp(-M^)  - /tTm  erfc  M},  the  exact  factor 
for  the  neutral  ion  flux.  Equating  Eqs.  (A-16)  and  (A-17),  we  may  write: 

-<j>  5 + ln2M^  + InvOnT  (A-18) 

1 e 

This  will  yield  an  overestimate,  since  the  ion  flux  is  actually  larger  than  that 
given  by  Eq.  (A-17).  For  M*4  and  m^/mg3l836  for  hydrogen,  we  obtain  23,  ver- 
sus the  self-consistent  value  (^--17  presented  in  Section  3.  For  M=10,  the  poten- 
tial is  <>— 109.  Hence,  for  large  M we  may  obtain  a good  estimate  by  keeping  only 
the  first  term  on  the  right-hand  side.  Then  the  potential  difference  generated 
by  the  flow  becomes  independent  of  the  temperature  and  may  be  estimated  by 

where  m^  is  the  ion  mass  and  v is  the  velocity  of  the  flow.  Thus,  the 
potential  difference  in  volts  may  be  estimated  by  0.00519*m(au)v^(km/sec)  where 
m(au)  is  the  ion  mass  in  atomic  units  and  v(km/sec)  is  the  flow  velocity  in 
km/sec.  Thus,  in  the  ionosphere,  assuming  m(au)»16  (O'*’  ions)  and  v»7  km/sec 
(orbital  velocity),  the  potential  difference  is  about  4.0  volts.  In  the  solar 
wind,  with  m(au)«l  (protons)  and  v»440  km/sec,  the  potential  difference  becomes 
1.0  kilovolt. 
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